Field Dependent Superfluid Density in the Optimally Doped 
SmFeAsOi-^F^ Superconductor 
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PACS 74 . 70 . Xa - Pnictides and chalcogenides 

PACS 74 . 25 . Ha - Magnetic properties including vortex structures and related phenomena 
PACS 76.75.+i - Muon spin rotation and relaxation 
PACS 74 . 25 . Bt - Thermodynamic properties 

Abstract. - The magnetic field dependence of the in-plane magnetic penetration depth A a (, for 
optimally doped SmFeAsOi-^Fy was investigated by combining torque magnetometry, SQUID 
magnetometry, and muon-spin rotation. The results obtained from these techniques show all a 
pronounced decrease of the superfluid density p s tx A~ 6 2 as the field is increased to 1.4 T. This 
behaviour of p s is analysed within a two-band model with self-consistently derived coupled gaps 
and p s = psi + p S 2, where p s i related to the larger gap is field independent, and p S 2 related to the 
smaller gap is strongly suppressed with increasing field. 



After the discovery of superconductivity in 
LaFeAsOi-xFz [I] with a transition temperature 
T c ~ 26 K, a whole new family of iron-based super- 
conductors was found with a maximum T c ~ 55 K for 
SmFcAsOi-^Fj, [2], where y < x due to a possible oxygen 
deficiency in the chemical composition [3J. Various 
experiments indicate multi-gap superconductivity within 
the family RE¥eAsO\- x ~F y (RE = rare-earth element) 
[3 8 . Furthermore, the magnetic penetration depth 
anisotropy, 7a = A c /A a ;,, increases with decreasing tem- 
perature, in contrast to the upper critical field anisotropy, 
j H — H^/Hjj^, which decreases with decreasing 
temperature [9HTT]. similar to those of the two-gap 
superconductor MgB2 [T21H3], although with reversed 
slopes [T4Tll6| . Here A, and denote the magnetic 
penetration depth and the upper critical field components 
along the crystallographic direction i (afo-plane or c-axis). 

Besides of the influence of the magnetic field H and 
temperature T on the anisotropy, the direct influence of 
H and T on A is essential in probing multi-gap supercon- 
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ductivity [T7]. However, it is difficult to obtain reliable 
experimental evidence for multi-gap superconductivity 
from the temperature dependence of A in samples contain- 
ing magnetic ions such as SmFeAsOi_ x F y . Importantly, 
the superfluid density p s = n s /m*, with n s being the 
superfluid carrier density and to* the effective carrier 
mass, can be probed directly by measuring A~ 2 cx p s . 

In the two-gap superconductor MgB 2 [T7] p s was 
found to be field dependent due to a suppression of 
the superfluid density in the band with the smaller gap 
[IH1IIS]. A similar dependence was observed for NbSc 2 
[mCH], V 3 Si YNi 2 B 2 C US], and La 1 . 83 Sro.i7Cu0 4 
[2~4] and again suggested to stem from either multi-gap 
superconductivity or to be related to changes of the 
internal field distribution of the vortex lattice. 

Previous studies of X a b in LaFeAsOi-xF^ [25] and in 
SmFcAsOi-a, and NdFeAsOi-a, [26] by muon-spin rota- 
tion (/xSR) experiments indicated a field dependent p s up 
to 0.6 T in the -R-EFeAsOi-zF^ family. Also in the related 
Ba(Feo.926Coo.o74)2As 2 compound a field dependent p s 
was reported for fields up to 0.2 T [27], and being inter- 
preted in terms of multi-gap superconductivity. However, 
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Fig. 1: Measured magnetic moment m of SmFeAsOi-^Fy as a function of temperature T in 1 mT magnetic field: a) single-crystal 
sample A, b) powder sample B. The rather sharp transition at T c demonstrates the good quality of our samples. 



high field investigations of SmFeAsOo.gFo.2 single crystals 
between 3 T and 30 T infer a field independent p s for 
SmPeAsOi-aFj, [28]. This obvious controversy is clarified 
in this letter by extended investigations of the magnetic 
field dependence of A a & in the SmFeAsOi_ x F y system. 
As the magnetic penetration depth X a b derived from ^tSR 
data must be carefully evaluated in order to demonstrate 
that the observed magnetic field dependence of p s is in- 
trinsic and not influenced by artificial contributions from 
the magnetic field distribution, additional experimental 
methods to investigate A a b are desired. Magnetisation as 
a global probe is less sensitive to the microscopic details 
of the magnetic field distribution inside a superconductor 
which is measured in pSR experiments. Additionally in 
magnetic torque experiments A a b is measured in terms of 
the average torque amplitude and almost insensitive to 
vortex pinning effects. Thus, a combination of magnetic 
and fiSK experiments is eligible to provide convincing 
experimental evidence for the magnetic field depen- 
dence of A a t>. Here, we report on a combined magnetic 
torque, SQUID magnetisation, and pSR study of the 
magnetic penetration depth, evidencing a pronounced 
magnetic field dependence of X a b in optimally doped 
SmFeAsOi-^Fy up to 1.4 T. 

In this work single crystal and powder SmFcAsOi-^Fj, 
samples were used, both of them synthesised by high- 
pressure techniques [3,29 . During the growth process, 
the composition SmFeAsOi_ x F y was fine tuned by 
varying x and y, in order to achieve a doping level close 
to optimal [3J. The single crystal of nominal composition 
SmFeAsOo.8Fo.2 (sample A) with an onset of T c ~ 47 K 
had approximate dimensions of 90 x 80 x 5 /im 3 . The 
powder sample of nominal composition SmFeAsO0.s5F0.15 
(sample B) with an onset of T c ~ 54 K was pressed 
into a pellet of cylindrical shape with 4 mm diameter 
and 2 mm height. Both samples were characterised by 
X-ray diffraction with similar properties as reported 
earlier [21 122] ■ Figure [1] displays the low- field magnetic 
moment m of sample A and B, measured in a SQUID 
magnetometer MPMS-XL (Quantum Design) using the 



reciprocating sample option (RSO) with pqH ~ 1 mT in 
zero-field cooled (zfc) mode [3j[29]. Both samples exhibit 
ideal diamagnetism at low temperature. The transition 
width (defined by the temperature range between 10% 
and 90% of ideal diamagnetism observed in the zfc 
low-field m(T) curve) is found to be 2 K and 6 K for the 
samples A and B, respectively, which demonstrates the 
good quality of the samples studied here. 

The single crystal sample A was investigated by 
torque magnctomctry in order to determine X a b(H). 
All measurements were performed using a custom-made 
torque magnetometer, previously employed to determine 
the anisotropic properties of tiny single crystals of iron 
pnictides [5J[T0]. The magnetic torque r = Ho(fh x H), 
related to the magnetic moment m, was recorded as a 
function of the angle 9 between the magnetic field H 
and the c-axis of the crystal. During the measurement, 
H was rotated clock- and counter clockwise around the 
sample at constant temperature T. Due to the present 
irreversibility, a vortex-shaking technique [9j I30J was 
applied with a 200 Hz electromagnetic field of a few mT 
perpendicular to H. The clockwise (6 + ) and counter- 
clockwise (t9~) torque scans were averaged according 
to T rcv (6») = [t(9+) + t(0")]/2, in order to further 
diminish irreversibility effects. A small sinusoidal normal 
state background of anisotropic paramagnetic origin was 
subtracted from all data by a symmetrising procedure 

[MIS]. 

A theoretical expression for the angular dependence 
of the torque, involving the two anisotropy param- 
eters 7a and 7# independently, was derived within 
Ginzburg-Landau theory by Kogan }31j : 
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Fig. 2: (Colour online) Reversible torque for the single crystal 
SmFeAsOi-zFy (sample A), measured at 36 K and 20 K in 
various magnetic fields. Solid lines are calculated using Eq. |T]), 
from which X a b was extracted. The inserts show the magnetic 
field dependence of \ ab . 



Here <!>o is the magnetic flux quantum and V the volume 
of the crystal. The scaling function €j(0) = [cos 2 (6>) + 
7~ 2 sin 2 (6')] 1 / 2 with i = A, H is different for 7a and jh, 
and e£(0) denotes the derivative with respect to 9. The 
torque data were analysed with Eq. (JlJ . From earlier work 
|9I10) the temperature dependence of was determined, 
and the anisotropy parameters were found to be field in- 
dependent up to 1.4 T. Since Eq. (TTJ) contains multiple 
correlated parameters, making a simultaneous analysis of 
all quantities difficult, all parameters except X a i, were kept 
fixed in the studied field range at the same values reported 
in [TO], thereby reducing the scattering of A a f,. Because, 
upper critical fields in optimally doped SmFeAsOi-^F^ 

are reported to exceed tens of Tesla (SHU [28], this re- 
lic 

striction to fix H" 2 for different H at the same tempera- 
ture during the data analysis appears reasonable since the 
field distribution of the vortex lattice at such low mag- 
netic fields is close to the London limit. Representative 
experimental torque data collected at two temperatures 
and various magnetic fields together with the calculated 
curves [Eq. |(TJ] are presented in Fig. [21 At 20 K and 36 K 
A a & increases substantially with increasing magnetic field 
(see inserts to Fig. [2J). For T < 20 K the presence of dis- 
turbing pinning effects [9] made the analysis meaningless, 
while for T > 36 K the torque signal was too weak to ex- 
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Fig. 3: Results of the SQUID investigation for sample B (pow- 
der) of SmFeAsOi-zFy. Upper panel: Field dependence of the 
magnetic moment m at 45 K measured for increasing and de- 
creasing magnetic field. Lower panel: Field dependence of A~ b 2 
derived by Eq. 



tract reliable values of A a f>, especially in the relevant low 
field regime. 

Additional SQUID magnetisation and /iSR investiga- 
tions were performed on the larger powder sample (sam- 
ple B). For anisotropic layered superconductors the effec- 
tive magnetic penetration depth A e ff of unoriented powder 
samples can be related to X a b via A e g — 1.31A a b |32j . For 
optimally doped SmFcAsOi_. T F.y we find values of 7a, ex- 
ceeding 7 at T c and increasing to 19 at zero temperature 
9 10], thus confirming this approximation and being in ac- 
cordance with other experiments which reveal comparable 
values for 7a in related compounds [55J[53]. Field depen- 
dent SQUID measurements of the magnetic moment at 45 
K are shown in the upper panel of Fig. [3] The analysis of 
these data is based on the average of the magnetisation for 
increasing and decreasing magnetic field with A~ b 2 being 
derived as [T7| 
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From Eq. ([2j! A^ 2 was determined and is shown in the 
lower panel of Fig. [3] Obviously, A~ b 2 at 45 K decreases 
continuously as the field is increased up to 1.4 T, in good 
agreement with the torque result (see Fig. 0]). 

In order to substantiate our main result on the mag- 
netic field dependence of A a t obtained on single crystals, 
the powder sample B was further studied by transverse- 
field (TF) and longitudinal-field (LF) /xSR at the Paul 
Scherrcr Institute (PSI) (For detailed information on the 
TF and LF /xSR techniques see e.g. [M]). Previous /xSR in- 
vestigations of i?SFeAsOi_ a; F y [26l[35H37] indicate that: 
(i) The muons stop at two lattice sites, one is close to 
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the magnetic moments of the Fe ions in the FeAs layers, 
whereas the other is near the spacing layers containing the 
rare-earth ions [26j[35j|36] . (ii) The magnetic moments of 
the Sm ions order antifcrromagnctically at temperatures 
below ~ 5 K, leading to a pronounced local-field broad- 
ening of the magnetic field distribution - In addi- 
tion, for temperatures above 5 K the finite correlation of 
the magnetic moments leads to thermally activated spin- 
fluctuations and causes a damping of the TF time spectra 
2li 137"]. (iii) Consistent with the above two points, the 
fiSR time spectra of SmFeAsOi_ :E F. y consist of a slow and 
a fast relaxing component J351I3B| . Combining these three 
points and invoking, according to calculations by Hayano 
et al. |38] , the magnetic correlation of the Sm moments to 
be similar in LF and TF experiments, the total transverse- 
field asymmetry A(t) of the time evolution of the muon- 
spin polarisation in SmFeAsOi-^Fy can be written as |26] 

A(t) = (A s e~ Ast + ,4 f c- Aft ) e^ 2 * 2 / 2 cos(s„Bi + 0). (3) 

Here ^4 S (^4f) and A s (Af) are the asymmetry and the de- 
polarisation rate of the slow (fast) relaxing component, re- 
spectively, B is the magnetic induction at the muon sites, 

is the muon gyromagnetic ratio, 4> is the initial phase of 
the muon-spin ensemble, and a is the Gaussian measure 
of the second moment of the magnetic field distribution 
within the sample. 

All /zSR time spectra were fitted simultaneously by 
Eq. §5§ with the ratio A s /A{ as a common parameter and 
A s and Af as inidividual parameters for each data set. The 
total asymmetry A s + Af was kept constant within each 
set of the data. The temperature dependent depolarisa- 
tion rates A s and Af in Eq. ^ are both assumed to have 
the form 



Ai(T) Ai(0) 



■Cexpi-Eo/keT), 



s,f, 



(4) 



with Eq being the activation energy for both relax- 
ations and C a constant. Our analysis yields A s (0) ~ 
0.20(5) /us" 1 , A f (0) ~ 2.71(7) /xs" 1 , and E ~ 20(1) meV. 
This procedure includes the background contribution due 
to the correlation of the Sm magnetic moments. Since 
the Sm moments order antifcrromagnctically below ~ 5 K 
26,315], an analysis for T ^ 10 K was omitted in this 
work. From Eq. ([3]) the two-component Gaussian relax- 
ation rate a was obtained: a = (of c + cr^m) 1 ^ 2 , where 
(T nm is the temperature independent nuclear moment con- 
tribution, and <7 SC is the superconducting relaxation rate 
due to the muon-spin depolarisation in the vortex lattice 
[25] . The values of X a b — A e ff/1.31 were derived from ct sc 
according to [39] 



<r sc [MS _1 ] = 4.83 x 10 4 (1 - H/H c2 ) 



(5) 
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with the same values of H C 2 as used for the torque ex- 
periment. The magnetic field dependence of A~ & 2 , derived 
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Fig. 4: (Colour online) Field dependence of A~ b 2 oc p s of 
SmFeAsOi-ajFy derived from the magnetic torque (sample A, 
closed symbols) and /xSR (sample B, open symbols) experi- 
ments. 



from the /iSR data together with the results from extended 
torque studies, are presented in Fig. |U Evidently, A~ b 2 is 
substantially suppressed with increasing magnetic field at 
all temperatures studied and follows, when normalised to 
Alt , (1.4 T/fio), on a universal curve (see Fig. [5]). 

The present results as obtained from all three tech- 
niques demonstrate very clearly a decrease of A~ b 2 with 
increasing magnetic field. From Figs. [3] [4] and[5]wc esti- 
mate a suppression of the superffuid density by 20% when 
the field increases from T up to 1.4 T. Assuming a single- 
gap scenario, for H <C H}}£ , only a nodal superconductor 
is expected to show a field dependent superfluid density 
[2T) ] l2% l l3"T p tJ]. Since however SmFeAsOi-^F^ is nodeless 
8,41], alternative explanations need to be considered. An 
influence of H c \ on the determination of can be ex- 
cluded, since for an estimated value of A a h(0) ~ 200 nm the 
corresponding H^(0) oc A~ h 2 (0) is of the order of 0.05 T, 
well below the lowest field studied here. A systematic 
error in the determination of A a b caused by pinning ap- 
pears to be unlikely as well: (i) In SQUID magnetometry 
pinning contributes strongly at low temperatures. The 
field dependence of the magnetic moment was performed 
at 45 K, where pinning effects are considerably reduced, 
(ii) In the /iSR experiments all data were obtained in field 
cooled vortex-lattice states. From the field dependence of 
the magnetic field distribution in the mixed state, a distor- 
tion of the fiSR spectra caused by pinning can be excluded 
[26j . (iii) In the torque experiment, a shaking technique 
was applied in order to minimise irreversibility effects |30] . 
Thus, the determination of A a f, is not affected by pinning. 
Furthermore, the derived X a b is mainly susceptible to the 
average torque amplitude which is almost insensitive to 
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Fig. 5: (Colour online) Universal behaviour of p s observed for 
both samples (A and B) of SmFeAsOi-^Fy for all investigated 
temperatures, normalised to its value at the maximum field of 
1.4 T. The black line is calculated using Eq. ((6J. The meaning 
of the symbols is the same as in Fig. U 



pinning, in contrast to H'J^, which was fixed during the 
data analysis [WTU]. 

As already mentioned above, the observed magnetic 
field dependence of A~ b is absent in nodeless single-gap su- 
perconductors [42] but has been observed in various multi- 
band superconductors, where the cases of MgB 2 fT7HT§] 
and Lai.83Sro.i7Cu04 [21] have been analysed in detail. 
The analogy to these systems suggests that related physics 
apply here as well. Assuming that SmFeAsOi-^F^ is a 
fully-gapped superconductor (HHT], a two-gap model [43] 
is used to calculate the superfluid density p s = p s \ + p S 2, 
stemming from the coupled gaps (Ai,A 2 ). In this ap- 
proach, the zero-temperature gaps are calculated self- 
consistently in relation to the total superfluid density 
p s . The best agreement is found by fixing the values 
Ai(0) = 13.83 meV and A 2 (0) = 5.26 meV, close to 
the s-wave gap values as obtained by point-contact spec- 
troscopy {6171, performed on samples from the same source 
[25] . With increasing magnetic field the superfluid den- 
sity is mainly suppressed in the band with the small gap 
(p s2 ,A 2 ) where the corresponding intra-band interaction 
approaches zero with increasing field. Note however, that 
the gap itself is not zero due to finite inter-band inter- 
actions analogous to MgB 2 [H1[T9]. The total superfluid 
density p s (H) is expressed as a sum of a contribution from 
the band with the large gap, p s \, which is independent of 
the field, and a field dependent contribution from the band 
with the small gap, p S 2 oc \j\fii. The field independence 
of psi stems from the same physics as in MgB 2 [T8lfT9] and 
the l/y/H dependence of p S 2(H) is based on the multi- 
gap analysis of Lai.ssSro 17CUO4, where the coupled gaps 
have different field dependencies, namely, the d-wave gap 
being proportional to \H [44] and the s-wave gap hav- 
ing a l/\/H dependence [45] . Since the s-wave symmetry 
of the SmFeAsOi-zFy superconductor is well justified [5], 
the 1/VH dependence is plausible. Thus, the expression 



p s (H)/p s i is derived as 

p s (H)/p sl = 1 + C(H/H )- 1/2 /psi, (6) 

where C denotes the amplitude of p s2 (i?) and Hq is the 
characteristic field for suppressing p s2 . A fit to the nor- 
malised superfluid density data by Eq. (jSJ) leads to an esti- 
mation of (C I Psi) 2 PqHq ~ 6.1 mT. The calculated result 
is shown in Fig. [5] in comparison with the normalised data, 
with obviously good agreement between both. 

It is important to emphasise that although p s is sup- 
pressed by 20% in small fields, T c is almost not changing 
on the same field scale, consistent with high values of H C 2. 
This can be understood in terms of a suppression of the 
intra-band coupling within the band of the small gap, hav- 
ing a pronounced influence on the corresponding p s , but 
almost no effect on T c [TMP^] . 

In conclusion, we observe a magnetic field dependence 
of the superfluid density in the iron-pnictide supercon- 
ductor SmFeAsOi-^Fj,. This experimental observation 
can be consistently described within a coupled two-gap 
model, where the superfluid density related to the large 
gap is field independent, whereas the one of the small gap 
is strongly suppressed with increasing magnetic field. To- 
gether with the temperature dependence of the anisotropy 
parameters, multi-gap superconductivity is well reflected 
in the thermodynamic behaviour of SmFcAsOi_ x F y . 
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